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Abstract 

The formation mechanism of the barium stars is thought to be well understood. Barium-rich material, lost in a stellar 
wind from a thermally-pulsing asymptotic-giant branch star in a binary system, is accreted by its companion main- 
sequence star. Now, many millions of years later, the primary is an unseen white dwarf and the secondary has itself 
evolved into a giant which displays absorption lines of barium in its spectrum and is what we call a barium star. A 
similar wind-accretion mechanism is also thought to form the low-metallicity CH and carbon-enhanced metal-poor 
stars. Qualitatively the picture seems clear but quantitatively it is decidedly murky: several key outstanding problems 
remain which challenge our basic understanding of binary-star physics. Barium stars with orbital periods less than 
about 4,000 days should - according to theory - be in circular orbits because of tidal dissipation, yet they are often 
observed to be eccentric. Only one barium-star period longer than 10 4 days has been published although such stars 
are predicted to exist in large numbers. In this paper we attempt to shed light on these problems. First, we consider 
the impact of kicking the white dwarf at its birth, a notion which is supported by independent evidence from studies 
of globular clusters. Second, we increase the amount of orbital angular momentum loss during wind mass transfer, 
which shrinks barium-star binaries to the required period range. We conclude with a discussion of possible physical 
mechanisms and implications of a kick, such as the break up of wide barium-star binaries and the limits imposed on 
our models by observations. 
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1. Introduction 

The barium stars are population I red giants which show 
strong absorption lines of barium in their spectra. They 
make up about 1% of all G/K giants. They are not evolved 
enough to have produced the barium themselves, rather 
it was made in a companion star which has long since 
ceased nuclear burning. This is borne out by the obser- 
vation that all barium stars ar e in binary stellar systems 
(|McClure fc WoodswortbL Il990h . 

Barium is an s-process element made during the 
thermally-pulsing asymptotic giant branch (TPAGB) stel- 
lar evolutionary phase. In sufficiently close binaries contain- 
ing a TPAGB star, the primary becomes large enough that 
Roche-lobe overflow (RLOF) terminates the TPAGB and 
lessens the production of barium, preventing barium-star 
formation. At large separations the strong stellar wind of 
the TPAGB star leads to mass accretion on the second- 
ary and this is the favoured channel for barium star forma- 
tion. An intermediate regime may exist in which some wind 
mass -transfer occurs before RLOF begins (e.g. lHan et all 
Il995h . 

While the wind-accretion scenario is generally accep- 
ted, it does not explain the distribution of eccentricit- 
ies and periods of barium stars. Most barium stars with 
P < 500 days are in (near-)circular orbits as predicted by 
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tidal circularisation theory (e.g. IZahnl . 119771 Il989f) . Most 
of the rest of the barium stars have perio ds between 500 
and 1 4 days and eccentricities < e < 0.4 (jJorissen et all 
EMI). 

Population synthesis studies remain the best way to 
study the orbital characteristics of bariu m stars in a stat- 
istical manner. Models such as those of iPols et al.l ([2003) 
confirm that our understanding of barium-star formation 
is incomplete. While the models confirm that short-period 
systems are indeed (almost) circular, they predict circu- 
larisation for all systems with periods shorter than about 
4,000 days, which is not seen in the observations. They also 
predict a long-period tail of eccentric barium stars extend- 
ing to 10 5 days, also not seen in the observations. While 
binaries with periods greater than 10 4 days are difficult to 
detect, the fraction of barium stars with measured periods 
is large (35/37 for strong barium stars) so a population with 
undetected, longer-period orbits can be ruled out as an ex- 
planation of the discrepancy between models and observa- 
tions. The combination of large eccentricities and relatively 
short periods of the barium stars remains an unresolved 
problem. 

If tides are as efficient as predicted a mechanism must 
exist which generates eccentricity in barium star sys- 
tems at the end of the TPAGB phase of the primary. 
This is evident from the observed eccentricities of bin- 
ary post-AGB stars which show a distribution strikingly 
similar to the barium stars. Many sources of eccentri- 
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city have been i nvesti g ated s uch as enhanced mass-transfer 
at periastro n rtSokerl. l2000ft. interaction with a circum- 
binary disc (jFrankowski fc Jorissenl 120071) and eccentricity 
umping induced by a wind-R LOF hybrid mass transfer 
Bonacic Marinovic et all , 12008ft . 

Recently and independently of the study of barium stars 
several authors have suggested that white dwarfs are kicked 
when they are born, i.e. at the end of the TPAGB phase. 
This result is based on a comparison of the velocity disper- 
sion of young white dwarfs to old white dwarfs and main- 
sequence stars in globular clusters (jDavis et all |2008|) . The 
young white dwarfs have a higher velocity dispersion than 
expected by a few kins -1 and are cited as a p ossible mech- 
anism to prevent g lobular cluster collapse (jHevll . l2008t 
iFregeau et all I2009D . If white-dwarf kicks occur they must 
affect the orbital parameters of the barium stars as well 
as their low-metallicity equivalents, the carbon-enhanced 
metal-poor (CEMP) and CH stars. 

The cause of the white-dwarf kicks is unknown. 
Candidate mechani sms include asymmetr ic mass loss dur- 
ing the AGB (e.g. iFellhauer et all l2003h . magnetic fields 
(jSpruid . 11998ft and perhaps non-radial stellar pulsations. 
Dynamical interactions inside the globular cluster may also 
play a role but it is not clear why these should affect 
only young white dwarfs. Studies of planetary nebulae may 
be able to constr ain any connection be tween asymmet- 
ric mass lo s s (e.g . ISahai fe Traugeii 11998ft . magnetic fields 
(|Lee et all l2007f ) and a kick to the central, proto-white 
dwarf, star. However, it is often difficult to determine the 
velocity of the central star compared to the expanding plan- 
etary nebula because of interaction of the nebula with a 
non-uniform interstellar medium. 

In this paper we investigate the possibility that a kick to 
the newly-born white dwarf is the cause of the eccentricity 
in the barium stars. In Section [5] we describe our popula- 
tion synthesis models, in Section [3] we present our results 
and compare them to observed barium star periods and 
eccentricities, Section @] discusses successes and potential 
problems of our model before we conclude with Section \5\ 

2. Modelling Ba stars with white-dwarf kicks 

O ur binary population synthe sis model is based on that 
of lHurlev. Tout, fc Polsl l|2002ft with nucleosynthesis as de- 
scribed bv llzzard et all (|2004 120061 . 12009T ) . Binary systems 
hav e initial primary masses Mi dis tributed according to 
the iKroupa. Tout, fc Gilmord" (|l993ft initial mass function 
between 1.2 and 3Mq, a flat distribution in the mass ra- 
tio q = M 1 /M 2 (where Mi > M 2 ) and a flat distribu- 
tion in log a between 400 and 10 5 Rq. Initial eccentricities 
are chosen from a distribution /(e) oc e between and 1. 
We assume that all stars form in binaries and our mod- 
els here have a grid resolution Nmi x _ZVm2 x N a x N e — 
30 x 30 x 60 x 30. 

We model tides following iHurlev. Tout, fc Polsl (12002ft 
who base their model on the formalisms o f lZahnl (jl977l ) and 
iHutJ (|1981l ). The tidal circularisation ti mescale duri ng the 
TPAGB r circ - 10 4 (a/3i?) 8 years (e.g. ISokeri l2000h . This 
implies that tides rapidly circularise any binary containing 
an AGB star with a separation less than a few stellar radii. 

Systems which are close enough to enter Roche-lobe 
overflow during the AGB undergo common envelope evol- 
ution with rapid circularisation, ejection of the stellar en- 
velope and little barium production and accretion onto the 



secondary. These systems lead to either little barium pro- 
duction or perhaps short-period barium stars. Some authors 
have speculated about stable case-C mass transfer which 
may lead to the formation of the shortest-period barium 
stars (e.g. lHan et ail . 11995ft and the intrig uing possibility 
of a radiation-distorted Roche geometry (jDermine et all 
2009). We do not explore this channel in detail here be- 
cause we focus on the wind mass-transfer scenario. 

Our nucleosynthesis mo del includes th i rd dr edge up 
with an efficiency given by iKarakas et all (12002ft and s- 
proces s abundances based on the models of iBusso et al.l 
(|200ll ). Our initial metalli city is Z = 0.008 wit h an a bund- 
ance mixture according to lAnders fc Grevessel (|1989| ) . The 
13 C pocket efficiency, 13 £, is a multiplicative factor used to 
change the amount of of 13 C relative to the 2.8 x 10 -6 M Q 
in the standard pocket of iGallino et al.l (jl998ft . We set it 
to 1 in order to make sufficient barium stars at Z — 0.008 
(about 1% of GK giants). At Z = 0.02 even with 13 £ = 2 
we cannot make enough Ba stars, however the amount of 
third dredge up is also rather uncertain and can be in- 
cr eased artificially t o obtain the desired 1% ratio (as noted 
bv lHan et all 119951 ). We do not try to constrain either 13 £ 
or the amount of third dredge up, this has bee n attempted 
by others (e.g. iBonacic Marinovic et all . 120061 see also the 
discussion in Section |4~4]) . 

Mass l oss during the TPAGB is p arameterised by the 
formula of IVassiliadis fc Wood! l|l993ft and mass accretion 
onto t he secondary is at the rate given by Bondi fc Hovlel 
as described in IHurlev. Tout, fc Polsl (|2002l Eg. 6) 



with an efficiency parameter aw = 1-5 and an accretion 



rate limited by 



Mi 



< 0.8 



Mi 



(where star 1 is the donor, 



star 2 the accretor). At the end of the TPAGB phase of the 
primary star a kick is applied to the nascent white dwarf 
with a fixed speed ctwd a nd random direction a c cordi ng to 
the method described in IHurlev. Tout, fc Polsl (|2002l ap- 
pendix Al, but without the Maxwellian velocity distribu- 
tion). We apply the kick whether the TPAGB is terminated 
by a stellar wind or common-envelope ejection. 

We count stars as giants after they leave the main se- 
quence (as defined by IHurlev. Tout, fc Polsl 120021) and be- 
fore they become white dwarfs. We select stars with only 
G an d K spectral types c o rresp onding to 3620 < T e g/K < 
5150 (|Jaschek fc Jaschekl Il995f) . Of these, strong barium 
stars have [Ba/Fe] > 0.5 and mild barium stars have 
0.2 < [Ba/Fe] < 0.5. 



3. Simulated Ba-star populations 

In this section we present the results of our population 
syntheses with the systematic inclusion of progressively 
speculative physics as required to reproduce the observed 
e — log P distribution of the barium stars. 

We compa r e our models with the observations of 
iJorissen et al.l ([1998) but without the triple system 
BD +38° 118. For the strong barium stars the sample 
is almost complete (35 out of 37 have measured periods 
and eccentricities). Our weak barium s tar sample is aug- 
mented with observations of 56 Uma (Griffin, 2008) and 
the latest results of the HERMES survey (Dermine and 
Jorissen, private communication). Very few barium stars 
do not have a measured period and eccentricity and the 
missing stars do not affect any of the general conclusions 
we present in this paper. 
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Figure 1. Our canonical barium-star e — logP distribu- 
tion. Darker greyscale/colour indicates a greater number 
density of systems. The plus s ymbols + are the obser- 
vations of iJorissen et al.l (|1998f ) without the triple sys- 
tem BD +38° 118 , the diagonal crosses x are the latest 
HERMES observations (Dermine and Jo rissen, privat e 
communication) and the square □ is 56 Uma (|Griffinl . l2008h . 
Where points are joined by lines there are multiple possible 
orbital solutions. 



3.1. Our canonical model 

Our canonical model result, i.e. without white-dwarf kicks 
or otherwise modified physics, is sh own in Fig. [T] and well 
reproduces previous stud i es suc h as iPols et all ([2003) and 
iBonacic Marinovic et al.l l|2008h . Our simulated barium- 
star systems with periods of less than about 4,000 days are 
circular. At longer periods our modelled barium stars are 
eccentric and some have periods as large as 10 5 days. The 
observed barium stars do not match either property of our 
simulated systems, most being eccentric and having periods 
typically less than 10 4 days. 

An obvious solution to the long-period tail is to decrease 
the wind accretion efficiency through the aw parameter. 
However, this has the effect of reducing the number ratio 
of barium to G/K giants to significantly below 1% so is 
not justified by the observations. It may be that a more 
complicated dependence of «w on the orbital separation 
allows the formation of barium stars only at periods less 
than 10 4 days. However, without a suitable formalism it is 
difficult to test anything but an even more ad-hoc solution 
to this problem than simply altering the constant a\v- As 
such, we keep aw = 1-5 in the remainder of this paper 
(except where stated otherwise). 

The gap in our simulated population at ~ 1,000 days is 
because of our chosen common envelope evolution prescrip- 
tion. The AGB core and companion star spiral together in- 
side the common envelope so the period shortens, in some 
cases to less than 100 days. Some observed barium stars in- 
deed have such short periods, but are not necessarily circu- 
lar. We discuss our common envelope prescription in more 
detail in Section I3~3l below. 
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Figure 2. As Fig.Q]but for a population including 4kms~ 1 
natal white-dwarf kicks. 
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Figure 3. As Fig. [3] but for a population with efficient 
common-envelope ejection (Ai on = 0.1). 

3.2. White-dwarf kicks 

In Fig. [2] we show our simulated barium-star population 
with the inclusion of a 4 km s _1 kick at the moment of birth 
of the white dwarf. Our simulation and the observations are 
in much better agreement than in the model of the previous 
section (c.f. Fig.[T]). Even the short-period systems acquire 
a small eccentricity. However, the problematic long-period 
tail extends to even longer periods than in our canonical 
model. These systems are weakly bound and a kick forces 
the binary into a wider orbit. However, too strong a kick 
leads to many broken binaries which are not observed, see 
Section O 



3.3. Common envelope evolution 

Our common-envelop e prescription is bas ed on the energy- 
balance formalism of iHurlev et al.l 1 20021) with a = 1 an d 
A fitted to detailed stellar models ( Dewi fe Taurisl . [2000) . 
It leads to the spiral in of systems with periods around 
1,000 days an d hence the perio d gap observed in Fig. [T] (as 
noted by e.g. lHan et ail Il995f ). The period gap is not seen 
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Figure 4. As Fig. [3] but for a population ver y efficient or- 
bital a ngular momentum loss according to the lHurlev et al] 
( 2002) formalism (Eq.ffland an accretion parameter aw = 
100, see text for details). 



in the observed barium-star period distribution and so one 
is forced to ask how it can be removed. 

The most obvious solution is to increase the efficiency 
of common-envelope ejection and hence reduce the amount 
of spiral- in before the envelope is lost. We do this by adding 
a fraction Ai on of the recombination energy of the envelope 
to the energy used to eject the envelope. In the simulations 
described in Sections 13.11 and 13.21 A; on = 0. However, be- 
cause TPAGB envelopes are cool and weakly bound, even 
a small value of Ai on , just a few per cent, is sufficient to 
close the period gap, as shown in Fig.[3]for Ai on = 0.1. For 
the shortest-period barium stars a smaller value of Ai on is 
required which suggests it may not take a constant value 
as we have assumed. W e note that use o f the common- 
envelope prescription of iNelemans fc Toutl (|2005| ) gives a 
similar result. 



3.4. Orbital angular momentum 

Wind accretion is not likely to be 100% efficient, that is, 
some material is always lost from the system. The material 
that is lost carries away orbital angular momentum so the 
orbital separation changes. It is not clear how much angu- 
lar momentum is lost per unit mass ejected. Our standard 
model, as presented in the se ctions above, uses the prescrip- 
tion of iHurlev et al.l (|2002f ) which provides the following 
formula, 



J, 



J, 



orb 



orb 



Mi + M 2 



Mi 



Ah 
Mi 



M 5 



(1) 



which is always negative because Mi < and M 2 > 
(where star 1 is the donor TPAGB star and star 2 is the 
accretor). When wind accretion is very efficient, the term 
— M2M1 is significant, the or bit shrinks and we obtain the 
results shown in Fig.SJ In the lHurlev et al] (|2002ft model ef- 
ficient wind accretion can be simulated with a large value of 

aw- The accretion rate is then limited to M 2 < 0.8 Mi 
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Figure 5. As Fig.0]but with our alternative orbital angular 
momentum loss formula (Eq. [5] with I = 2 and an accretion 
parameter aw = 1-5). 



Figure [1] shows the case with aw = 100 which is effect- 
ively 80% efficient accretion because of the imposed limit 



M 2 



< 0.8 



Mi 



Most of the long-period barium stars in 

our model shrink their orbits to periods less than 10 days 

so agree much better with the observations . 

One may ask whether the IHurlev et al.l (|2002f ) formal- 
ism should be believed. If Mi = —M 2 , i.e. mass transfer is 
conservative, one would expect J or b = which is clearly not 
the case. This has lead us to consider a simpler prescription 
in which material lost from the binary system carries away 
some multiple of the specific orbital angular momentum, 



(2) 



where I is a free parameter. This formula has the advantage 
that J or b = in the case of conservative mass transfer. It is 
difficult to constrain I but Fig. \5\ shows the case I — 2 which 
looks very similar to Fig. 0] but without the requirement of 
a very large aw- The long-period tail is always present, but 
the number of stars we predict there is small. 



4. Discussion 

The physical mechanism for the white dwarf kic k is un- 
know n but candidates includ e magnetic fields (|SpruilJ . 
Il998f ). asymmetric mass lo ss (iFellhauer et all l2003h and 
mass transfer at periastron (|Sokerl . l2000l ). According to the 
evidence from globular clusters the m echanism must be at 
work in both single and binary stars (j Davis et all . 120081) . 

We are working on an alternative mechanism which may 
be of relevance to barium stars, that of eccentricity pump- 
ing due to a circumbinary disk (Dermine et al. in prep.), but 
again this cannot be active in single stars. It would seem 
that asymmetric AGB wind loss and/or interaction with a 
magnetic field is the leading candidate if white-dwarf kicks 
are to occur in both single and binary stars. Still, a number 
of open questions remain and we briefly address a few of 
them here. 
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Figure 6. The fraction of binary systems which are disrup- 
ted as a result of a white-dwarf kick as a function of the 
kick speed. The models are as those of Section[32](Ai on = 0, 
aw = 1.5 and orbital angular momentum loss according to 
Eq.Q). 



4.1. Fast and slow kicks - and in which direction? 

In all our models presented above we have assumed an in- 
stantaneous kick is given to the white dwarf at the end of 
the primary AGB phase. However, any kick which is impar- 
ted by an asymmetric wind is likely to be effective during 
the superwind phase of the primary AGB star, i.e. on a 
timescale of ~ 10 4 years. This is of the same order as the 
tidal circularisation timescale and much longer than the 
orbital period. A coupling of the orbital dynamics with a 
"slow" kick is beyond the scope of this paper but there may 
still be an effect on the eccentricity particularly if the kick 
is coupled to the orbital phase. 

The direction of the kick is also an issue crucial to the 
further evolution of the binary system. Close barium stars, 
with P < 4,000 days, are circular by the end of the AGB 
so any kick increases the eccentricity of the orbit. In wide 
systems which remain eccentric throughout the AGB phase 
the eccentricity may increase or decrease depending on the 
kick direction. For simplicity, we choose a fixed kick speed 
and a random direction: this may not be realistic in a bin- 
ary system in which there is a clear axis of symmetry. The 
presence or otherwise of a preferred direction would depend 
on the kick mechanism, which is still unclear. 

4.2. Disrupted systems 

A consequence of kicking a component of a binary star 
which is well known in the study of massive stars is disrup- 
tion of the binary system. The fraction of binaries which 
are broken as a result of kicks in our otherwise canonical 
model (c.f. Section [5]) is shown in Fig. H2 For small velo- 
cities, < 3kms , the number of broken systems is less 
than a few per cent . This is a useful con straint on the kick 
mechanism because [jorissen et al.l (|1998() present 65 binary 
barium stars, and possibly one single star in a complete 
strong barium star sample, hence the disruption fraction 
is less than one in 65, or about 1.5%, corresponding to 



a maximum ctwd = 2 — 3kms _1 . To the authors' know- 
ledge no conclusively-single giant barium stars have yet 
been found, but their number may be few and they may 
have been missed. Some may be present among subgiant 
CH stars. 

Wide binaries are most easily disrupted for a given 
cwd- We find a smaller breakup fraction when we con- 
sider our models with efficient orbital angular momentum 
loss because these barium stars have shorter periods and 
so are more tightly bound. For example, using Eq. [5] with 
I = 4 we obtain a breakup fraction of only 0.7% with 
cwd = 4kms _1 compared to 10% in our canonical model 
with the same kick speed. 

4.3. Mass loss, the angular momentum budget and Mira 

We have hitherto ignored the uncertainty introduced by 
our choice of mass loss rate. Our chosen rate, that of 
IVassiliadis fc Wood (1993), tends to be small during most 
the TPAGB but strong during the final superwind phase. 
This is advantageous for barium star formation as it allows 
for many thermal pulses and associated barium production. 
A smaller but more steady mass-loss rate may lead to fewer 
and weaker third dredge up episodes because the envelope 
is more quickly reduced in mass, but conversely any barium 
dredged up is less diluted so the amount of barium in trans- 
ferred material may be similar. In any case, the efficiency 
of accretion is so poorly understood that the abundance 
of barium in the primary is unlikely t o be the largest un- 
certa i nty in our analysis (a lthough see lStancliffe fc Jeffervl 
l2007t ICristallo et all [20091 for suitable analyses). 

It is not unreasonable to assume efficient loss of orbital 
angular momentum, e.g. Eq. [5] with I > 1, if there is signi- 
ficant locking of the stellar wind with a magnetic field or 
mass is lost from a circumstellar o r circumbinary disc . An 
alternative explanation is that the iHurlev et al.l (|2002h an- 
gular momentum formula, our Eq. [TJ is correct in the case 
of efficient accretion on the companion. Efficient accretion 
may already have been seen in systems such as Mira in 
whic h the stellar wind is c hannelled onto the companion 
star (Karovska et al., 2005). Recent simulations of wind- 
RLOF hybrid mass transfer suggest accre tion is indeed ef- 
ficient (|Mohamed fc PodsiadlowskH |2007[ ). An alternative 
int o which we are l ooking is the implementation of the work 
of ISepinskv et al.l (|2009h which would allow us to follow 
non-conservative mass transfer in detail, but this is beyond 
the scope of this paper. 

4.4. Mild vs strong barium stars 

The differences between the distribution of strong and mild 
barium stars in the e — log P diagram may help us constrain 
uncertain physics. Our models predict that strong barium 
stars are more likely to be made in the intermediate period 
range 10 3 — 10 4 days because at short periods RLOF trun- 
cates the primary AGB evolution while at long periods wind 
accretion is not efficient. In contrast, mild barium stars are 
made over the whole period range of 10 2 — 10 5 days. The 
iJorissen et all (| 19981 ) data may support this view. For the 
30 mild barium stars with known periods, the mean period 
is 2900 days with a standard deviation of 1900 days, while 
for the 35 strong barium stars with known periods the mean 
period is 2000 days with a standard deviation of 1700 days. 
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The new HERMES data further lengthens the mean period 
of the mild barium stars. 

There are five strong (binary) barium stars with peri- 
ods less than 500 days which prompt questions about their 
origin. It is possible that the presence of a companion star 
enhances third dredge up through the effects of tidal lock- 
ing and rotational mixing (the stellar core will not spin at 
the same rate as the envelope, we are working on detailed 
models of AGB binaries to test this hypothesis). Also, mass 
accretion may occur during the common envelope phase 
and l ead to strong barium star formation (jRicker 8z Taaml . 
l2008h . 

In principle we could use the properties and number of 
barium stars to constrain AGB physics, such as the amount 
of third dredge up, the efficiency of the 13 C pocket and the 
amount of s-process element production. In practice this is 
difficult because we can increase the amount of third dredge 
up in but decrease the 13 C pocket efficiency in our models 
yet obtain the same barium/GK giant fraction. Similarly, 
if we artificially increase the amount of dredge up in the 
primary AGB star we can reduce the efficiency of accretion 
on the secondary and obtain the same result. F or these reas- 
ons we have chosen the standard 13 C pocket of lBusso et al.l 
(l200lh with th e efficie ncy of third dredge up as predicted by 
iKarakas et~aTl (|2002D . At a metallicity Z = 0.008 we find 
that the Ba/GK-giant number ratio is ~ 1%, as observed. 

Simila rly, we could perhaps u se the masses of barium 
stars (e.g. lAllen fc BarbuvL l2006( ) to constrain the amount 
of accreted material, but without prior knowledge of the 
abundance of the accreted material, the amount of (or 
lack of) thermohaline mixing or the initial mass of the 
(now) bariu m star, this is diffi cult if not impossible, for 
exam ple see Husti et al.l (l2009f) and the CEMP star equi- 
valent iBisterzo et alJ (|200£ 



4.5. Implications for planetary nebulae 

Stars leaving the AGB pass through a planetary nebula 
phase before they go on to become white dwarfs. Any asym- 
metry in the AGB superwind might manifest itself as an 
asymmetric planetary nebula and indeed most planetary 
nebulae have an axisymmetric rather than spherical struc- 
ture. A large number present large-scale deviations from 
axisymmetry which can be explained by the presence of a 
binary companion. 



iSoker et al.l (|1998l ) investigated how an ec centric binary 
can explain deviations from axisymmetry and Soker (2000) 
proposed that the origin of the eccentricity is enhanced 
mass loss at periastron passages. However, this process re- 
quires a seed eccentricity which is not present in short- 
period barium-star progenitor systems which are circular- 
ised by tides. 

As white-dwarf kicks apparently occur in both single 
and binary stars, we should see the effect of the kick in 
single-star planetary nebulae. However, it is difficult to in- 
fer a kick from the structure of a planetary nebula, e.g. 
a shift in the position of the central star relative to the 
nebula, because its distortion can arise from the interac- 
tion between the nebula and the interstellar medium or the 
proper motion of the central star. 



5. Conclusions 

The combination of a white-dwarf kick with a speed of a 
few kms' 1 , efficient common-envelope ejection and orbital 
angular momentum loss enables us to much better match 
the distribution of eccentricities and periods of observed 
barium stars with our models. A small kick speed is compat- 
ible with observations of both globular clusters and barium 
stars. Efficient common-envelope ejection is required to pre- 
vent barium-star systems from spiralling to short periods 
during the common envelope phase. 

Our orbital angular momentum loss mechanism is less 
certain. With the exist i ng for malism in our binary model, 
based on iHurlev et alJ (|2002f ). we can shrink barium-star 
orbits into the observed period range with an accretion ef- 
ficiency of 80%. However, we are probably extrapolating 
their wind accretion model beyond its validity. If, instead, 
we assume that material lost from the system carries away 
twice the specific orbital angular momentum we achieve a 
similar result. 

In any case, a white-dwarf kick is a plausible mechan- 
ism for generating the eccentricity observed in the barium 
stars. It is supported by independent evidence from glob- 
ular clusters and is apparently not a uniquely binary-star 
phenomenon. If it acts in barium stars it should also be 
of consequence to the formation of CH and CEMP stars. 
Hopefully, future studies of asymmetric planetary nebu- 
lae, circumbinary disks in (post-)AGB stars, and globular- 
cluster and field white dwarfs may pin down both the origin 
of the white-dwarf kicks and barium-star eccentricity. While 
white-dwarf kicks may explain the barium star mystery, the 
evidence as it stands is circumstantial and concrete proof 
will only come from a combination of observations and im- 
proved modelling. 
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